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The global dynamics of gene regulatory networks are known to show robust- 
r*i ■ ness to perturbations in the form of intrinsic and extrinsic noise, as well as 

mutations of individual genes. One molecular mechanism underlying this ro- 

P^ ■ bustness has been identified as the action of so-called microRNAs that operate 

> '• 

t^. . via feedforward loops. We present results of a computational study, using the 

O ' 

O ' modeling framework of stochastic Boolean networks, which explores the role 

00 ! 
/ ■ that such network motifs play in stabilizing global dynamics. The paper intro- 

^ \ duces a new measure for the stability of stochastic networks. The results show 

that certain types of feedforward loops do indeed buffer the network against 

kJ ■ stochastic effects. 



The term canalization was coined by the geneticist C. Waddington to describe 
the theory that embryonal development is buffered against genetic and envi- 
ronmental perturbations. It is only recently that a molecular basis for this 
phenomenon has been suggested. Recent research has highlighted how the in- 
trinsic stochasticity of gene expression can drive changes in phenotypes. Short 
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segments of single-stranded RNA, so-called microRNAs (miRNA), represent an 
entirely novel agent of gene regulation discovered relatively recently, and have 
been proposed to function as canalizing agents that buffer the effects of such 
stochasticity in gene expression. According to this theory, when miRNA expres- 
sion is perturbed, stochasticity in gene expression can result in transitions to 
distinct cellular phenotypes. As miRNAs bind to gene targets they downreg- 
ulate translation of target mRNA into protein. Embedded in several different 
types of so-called feedforward loops (FFLs), miRNAs help smooth out noise and 
generate canalizing effects in gene regulation by overriding the effect of certain 
genes on others. 

Much experimental work remains to be done in elucidating this concept, and 
recent years have seen an explosive growth of publications in this area. There 
have also been a number of computational studies focused on canalization. In 
this paper, we carry out a computational study of the ability of the feedforward 
loop motif to buffer a gene regulatory network against intrinsic noise. This is 
done using stochastic Boolean network models as a computational instantiation 
of gene regulatory networks. We introduce a measure on networks that captures 
its "distance-to-deterministic" characteristics in terms of the stability of their 
attractors. For a given network, we successively introduce feedforward loops 
and track the resulting change in dynamics. The results show clearly that the 
feedforward loop motif buffers the network phenotype, in terms of stability of 
attractors, against perturbations from intrinsic noise. 

I. Introduction 

The term canalizaton was coined by the geneticist C. Waddingtoni to describe the theory 
that embryonal development is buffered against genetic and environmental perturbations. It 
is only recently that a molecular basis for this phenomenon has been suggested. Recent re- 
search has highlighted how the intrinsic stochasticity of gene expression can drive changes in 
phenotypes^. Short segments of single-stranded RNA, so-called microRNAs (miRNA), repre- 
sent an entirely novel agent of gene regulation discovered relatively recently^, and have been 



proposed to function as canalizing agents that buffer the effects of such stochasticity in gene 
expression 5,6 . According to this theory, when miRNA expression is perturbed, stochasticity 
in gene expression can result in transitions to distinct cellular phenotypes. As miRNAs bind 
to gene targets they downregulate translation of target mRNA into protein. Embedded in 
several different types of so-called feedforward loops (FFLs), miRNAs help smooth out noise 
and generate canalizing effects in gene regulation by overriding the effect of certain genes on 
others 7 -. Complex networks (viewed as graphs) ranging from the transcriptional networks in 
yeast and E. coli to engineered systems are enriched for certain graph-theoretic motifs that 
include feedforward loops^. 

An understanding of canalization in evolutionary biology is important as a cornerstone 
of natural selection and the emergence of new phenotypes 9 -, as well as for the understanding 
of diseases such as cancer. Transitions to new phenotypes have been implicated as one of 
the driving forces of tumorigenesis^~— ; and, interestingly, significantly altered expression 
of miRNAs is a feature of several cancers. Much experimental work remains to be done 
in elucidating this concept, and recent years have seen an explosive growth of publications 
in this area. There have also been a number of computational studies focused on canal- 
ization. Several papers have studied computational models that capture the evolution of 
canalization in networks and their ability to support significant mutation without change 
in the phenotype^ 1 ^, while others have studied models of stochastic gene expression as the 
internal source of noise in regulatory networks 17 . A detailed stochastic model of the abil- 
ity of miRNAs to buffer gene expression noise in a feedforward loop has been proposed 18 , 
providing evidence that this type of network motif can in fact play a canalizing role. There 
is evidence that miRNAs do not regulate their target genes directly; rather they act as 
post-transcriptional regulators by reducing the amount of mRNA and by repressing mRNA 
translation 19 , e.g., by binding to the 3'-UTR of a mRNA, which prevents this mRNA from 
being translated. 

In this paper, we carry out a computational study of the ability of the feedforward loop 
motif to buffer a gene regulatory network against intrinsic noise. This is done using stochas- 
tic Boolean network models as a computational instantiation of gene regulatory networks. 
We introduce a measure on networks that captures their "distance-to-deterministic" char- 
acteristics in terms of the stability of their attractors. For a given network, we successively 



introduce feedforward loops and track the resulting change in dynamics. The results show 
clearly that the feedforward loop motif buffers the network phenotype, in terms of stability 
of attractors, against perturbations from intrinsic noise. 

II. Modeling framework 

A. Stochastic Discrete Dynamical Systems 

In this study, the recently developed framework of stochastic discrete dynamical systems 
(SDDS) 20 is used to model gene regulatory networks. This framework is an appropriate 
set-up to model the effect of intrinsic noise on network dynamics. A stochastic discrete 
dynamical system in the variables xi, . . . , x n , which in this paper represent genes, is defined 
as a collection of n triplets 

where 

• fi : {0, l} n — > {0, 1} is the update function for Xi for all i = 1, . . . , n, 

• p] € (0, 1] is the activation propensity, 

• pi G (0, 1] is the degradation propensity. 

-s I 

The stochasticity originates from the propensity parameters p\ and p\, which should be 
interpreted as follows: If there would be an activation of Xj at the next time step, i.e., Xi(t) = 
0, and fi(xi(t), . . . ,x n (t)) = 1, then Xi(t + 1) = 1 with probability pj. The degradation 
probability pj is defined similarly. 

All variables are synchronously updated and one time step can be interpreted as the av- 
erage time needed for transcription and translation of the fastest of the genes considered. 
The propensity parameters for this fastest gene will be set to 1, and the propensity param- 
eters of genes with slower transcription and translation take proportionately lower values. 
Thus, this framework can be interpreted as introducing a very general stochastic sequential 
update scheme, which also allows for a variable not to be updated at all at a given step, a 
generalization of the usual approach in, e.g.,—. 



B. Quantifying Stochasticity 

In a deterministic discrete dynamical system, each initial configuration lies in exactly one 
basin of attraction. This changes when stochasticity is introduced. Now, from one initial 
configuration, different attractors may be reached. To measure the degree of stochastic- 
ity in particular dynamics, we look at every initial configuration and regard its transition 
probabilities to the different attractors. If each initial configuration only transitions to one 
attractor, the dynamics are deterministic, whereas lower maximal transition probabilities to 
attractors lead to proportionately more stochastic dynamics. In our context, the different 
attractors may be interpreted as different cellular phenotypes, which makes the connection 
to phenotypic stability discussed in the introduction^. Thus, this computational project fo- 
cuses on the stability of attractors under intrinsic noise, as it is affected by the introduction 
of feedforward loops. 

Based on this idea, we can define the degree of stochasticity in the dynamics of an SDDS 
F = (fiiPi,Pi)i=i- Let A(F) be the set of all attractors of F. Then we can calculate the 
average maximal transition probability to an attractor, where all 2 n state space configurations 
are considered and weighted equally, as follows: 

MF) = 1 £ (maxP(xA...AA))G[0,f]. (2) 

x€{0,l}" 

When F is a deterministic system, fi(F) = 1 always holds true. In comparison, for a 
stochastic system with a attractors, values as low as 1/a may be obtained; for stochastic sys- 
tems with a single attractor, fi(F) = 1 because the single attractor is eventually approached 
from any initial configuration. 

Most limit cycles that are attractors in a deterministic system are no longer attractors 
in an SDDS, because a cycle can be exited with a certain probability. Nevertheless, one 
particular kind of limit cycle, which consists of 2 fc elements and in which all but k bits are 
fixed, remains an attractor even in a SDDS. One such example is a 2-cycle, in which just one 
bit switches states, e.g., 000 -H- 001. Table [I] shows that such limit cycles occur very rarely 
by chance, and for computational reasons, we decided to include only steady states and limit 
cycles of length 2 5 or less in this study. This restriction does not influence the study since 
longer limit cycles that remain attractors in the SDDS practically do not occur. 



TABLE I. This table shows the average number of steady states and limit cycles that remain 
attractors in the SDDS framework for different network sizes. We found fewer than thirty such 
8-cycles among more than 250,000 networks of different sizes, and no 16-cycles at all. This shows 
that including only limit cycles of length 16 and less is not restricting the study. 



Network Size 


5 


15 


30 


50 


Cycle Length 


1 


2.8351 


3.6577 


4.3492 


4.8709 


2 


0.1529 


0.1522 


0.1540 


0.1486 


4 


0.0244 


0.0378 


0.0415 


0.0415 


8 


0.0002 





<0.0001 


<0.0001 


16 














Sample Size 


120000 


42000 


40000 


62500 



The basic procedure underlying the computational study is, for a given SDDS, referred 
to as the "basic" network, to construct several "extended" networks by successively adding 
nodes, which we shall refer to as miRNAs, together with one or more feedforward loops involv- 
ing the new miRNAs in a specific way. We then measure the change in the stochasticity mea- 
sure described above. Let F = (/j,pj,pt)rii be the basic network, and let F* = (/*, qj, qf)^ 
be the extended network, n\ < n<i. We hypothesize that the dynamics in the feedforward 
loop enriched network are less stochastic. To compare the dynamics of two systems with 
respect to their degree of determinism, we consider their difference m in /x-values 



m{F, F*) 







m(F«)- M (F) 



iif,(F*),fi(F) = l 
otherwise. 



(3) 



l-mm(fi(F),fj,(F*)) 

The denominator scales this difference into the range [— 1, 1]. Here, m = means that both 
networks have equally stochastic dynamics. If m is positive, the extended network F* is 
dynamically less stochastic than the basic network F, and a negative value of m suggests the 
opposite. The magnitude describes the difference in degree of stochasticity. A magnitude 
of 1 means that one of the networks has the dynamics of a deterministic system, whereas a 
magnitude of 0.5 suggests that one system is 50% less stochastic than the other. 



III. Methods 

For each set of input nodes we generated a certain number N of basic Boolean SDDS F, in- 
troduced miRNAs, in a way that will be specified later in this section, to obtain the extended 
network F*, and then compared their degree of determinism via m(F, F*). We considered 4 
network sizes n: 5, 15, 30, 50 nodes. The corresponding number iV is 20000, 7000, 5500, 5500, 
respectively. The networks generated are random, subject to the following restrictions. 

Large-scale studies of B. subtilis, E. coli and S. cerevisiae strongly suggest that the in- 
degrees of nodes in a transcriptional regulatory network are Poisson distributed with a mean 
of about two^. Thus, we chose a Poisson distribution with parameter 7 = 2 for the basic 
network. The only other restriction is that each gene is regulated by at least one other gene, 
which raises the average in-degree to approximately 2.2. The regulators of each gene are 
randomly chosen from the set of all n genes in the network, allowing self-regulation. 

All propensity parameters for transcription factors and miRNAs are also randomly chosen 
from [0.2, 1]. For computational reasons, the full interval [0, 1] is not used since a propensity 
parameter close to might strongly decelerate convergence to attractors, slowing the per- 
formance of the simulation. However, 0.2 as the lower limit for the propensity parameters 
still allows one process to happen up to five times more frequently than another. Each gene 
is regulated by a certain number of genes, depending on its in-degree, and random Boolean 
functions that actually depend on all input variables are used as update functions. 

After creating the basic network, we generate an extended network in a way reminiscent of 
post-transcriptional regulation by miRNAs. (Since we do not include a corresponding protein 
node for each gene node, the analogy is limited). Starting with the basic SDDS, miRNAs 
are iteratively added to this system by randomly choosing one gene as a transcription factor 
(TF) that induces transcription of the introduced miRNA. This miRNA, in turn, reduces the 
mRNA level of its own transcription factor; one example for such coregulation is the interplay 
between miR-133b and Pitx3 in midbrain dopamine neurons^. A lower TF mRNA level leads 
to a lower TF protein level, which then affects the regulation of all TF target genes. In the 
Boolean framework, a gene-specific threshold is used to distinguish between low (0) and high 
concentration (1). For some target genes, even after the TF mRNA reduction, there might 
still be enough transcription factor so that the target concentration is on the same side of the 



TABLE II. Example of how an additional miRNA is embedded as an input variable in a target 
gene's update function. The light grey rectangle contains the original random update function, 
in which 01, . . . ,04 € {0, 1} can be any Boolean values, with the sole restriction that the update 
function must depend on both inputs. Unless both miRNA and transcription factor are present, the 
miRNA has no influence since there is no TF mRNA that can be degraded or there is no miRNA 
that can catalyze this degradation. Only if both are present (gray rows), can the miRNA reduce 
the TF mRNA level to an extent that the target concentration changes because of the miRNA. 
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threshold as if no reduction had taken place; for other target genes, the target concentration 
might change significantly because of the TF mRNA reduction. Since this reduction is caused 
by the miRNA presence, the miRNA becomes a new regulator of the target genes in the latter 
case. One input variable in this study, called miRNA strength, describes the probability that 
transcription factor-target gene pairs fall into this latter case, i.e., that the TF mRNA level 
is significantly reduced, so that the Booleanized target concentration is the same as if no 
transcription factor had been present at all. If, for instance, the miRNA strength is 0.5, 
any miRNA regulates on average half of its transcription factor's target genes. However, we 
require each miRNA to regulate at least one target gene. This restriction ensures that each 
miRNA is part of at least one feedforward loop, consisting of transcription factor, miRNA, 
and target gene. Table [III depicts an example of how the update function of target genes 
regulated by a miRNA is expanded, taking into consideration the mode of action of miRNAs. 
Only if miRNA is present and TF mRNA is transcribed, the mRNA reduction takes place. 
In this case we see the same output as if no TF mRNA were present. 

Because all networks with exactly one attractor already have totally deterministic dy- 
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namics, in the sense denned earlier, i.e., \x = 1 for those networks, only basic networks with 
multiple attractors are considered in this study. Those are the interesting networks, in which 
actual stabilization of the dynamics might be observed. Particularly interesting dynamics 
occur if at least two attractors possess a relatively large basin of attraction. The network 
selection process therefore favors networks with multiple large basins of attraction by pick- 
ing a network only if at least two attractors are found more than once starting from twenty 
random initial configurations. 

If the extended network F* does not possess multiple attractors, m(F, F*) = 1 by defini- 
tion. One could argue that the loss of attractors in the extended network is one feature of 
stabilization through feedforward loops. On the other hand, however, this could be seen as 
an experimental bias. To consider both views, we use m to define two output measures, vfi\ 
and 7B2, one regarding any extended network and the other considering only those network 
pairs in which the extended network also possesses multiple attractors with at least two large 
basins of attraction. For a given set of input variables, we generate iV basic and extended 
networks, and measure 

N 

N 



1 v—^, 

mi = — y^{m(F, F*)|basic network F and extended network F* have multiple attractors}, 



(4) 

N 
1 - 

m 2 



— \^{yn{F, .F*) I only the basic network F is required to have multiple attractors}. 

(5) 



N 
i=i 



For any set of input variables, we expect m 2 > Vfi\ since all network pairs with less than two 
attractors in the extended network, which are omitted in m 1; have \x = 1 and thus a mean 
value closer to 1. However, we do not want to prefer one or the other measure and thus we 
report results for both, which have been obtained independently, i.e., a network pair that 
was used for m\ is not used for m 2 . 

A full analysis of the state space of a SDDS is only possible for small networks, so we 
used random sampling of initial configurations and an estimate of transition probabilities to 
attractors to approximate m(F,F*). We created a set of 100 random initial configurations, 
which were used in both networks to find the transition probabilities to attractors by up- 
dating each configuration 50 times, until an attractor was reached. In a small preliminary 
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study, we found that these two values yield a good trade-off between accuracy and efficiency. 

IV. Results 

Overall, we created over 300, 000 pairs of basic and feedforward loop enriched networks. 
The results for networks with sizes ranging from 5 to 50 genes can be seen in Figure [TJ None 
of these networks were required to be strongly connected, and in all of them the introduced 
miRNAs had full strength, meaning that each miRNA regulates all of its transcription factor's 
target genes in a feedforward loop structure. The main result is that both measures, m,\ and 
?772, are indeed positive for all network sizes and numbers of introduced miRNAs, indicating 
that miRNA-mediated FFLs can actually stabilize networks. It can also be seen that the 
impact of a single miRNA/FFL decreases when the network becomes larger. This means 
that larger networks require more miRNAs/FFLs for the same degree of stabilization. 

Table HTT1 shows the results for networks of size 50. We see that more miRNAs and thus 
more FFLs stabilize the dynamics. Whereas one miRNA of full strength with m\ ~ 4% only 
has a small impact, five such miRNAs already lead to mi ~ 12%, and the introduction of 
thirty miRNAs of full strength stabilizes the stochastic system quite a lot (mi ~ 0.37%). 
As expected, m% yields higher values and thirty miRNAs already reduce the stochasticity 
in the dynamics by more than 50% (772,2 ~ 0.57%). In the case that each miRNA regulates 
on average only half its transcription factor's target genes (but at least one), all values are 
considerably lower; the general behavior does not change, however. 

These results raise the question whether a given network can be fully stabilized by in- 
troducing a sufficient number of miRNAs. Indeed, under certain conditions this is possible 
by ensuring the existence of a unique steady state. If an n-gene network contains no self- 
regulating genes, then n miRNAs with full strength, each regulated by another gene, suffice 
to have fully deterministic network dynamics. Since the miRNA has full strength, it will 
downregulate any present mRNA, which ensures that only the value ai G {0, 1} (compare 
Table HV|) can be taken by the target gene at a steady state. Each gene is regulated by at 
least one other gene. Hence, each gene and its regulated miRNA can only take the value a\ 
in its truth table and the existence of a unique steady state is guaranteed. Thus, fi = 1 for 
such networks, which is equivalent to fully deterministic dynamics in the sense defined in 

10 



TABLE III. Comparison of the degree of stochasticity via mi and rri2 for not necessarily strongly 
connected networks of 50 genes, in which various numbers of miRNAs with full strength (Part 
a) and with strength 0.5 (Part b) are introduced. Overall, the more miRNA-mediated FFLs are 
introduced, the less stochastic the network dynamics become. 
a) 



Number 


of miRNAs 


1 3 5 8 10 15 30 


Average 


Number of FFLs 


2.401 7.298 12.36 19.79 24.77 37.27 74.50 


mi 


0.0374 0.0700 0.1083 0.1548 0.1756 0.2609 0.4042 


m 2 


0.0687 0.1444 0.2178 0.2832 0.3203 0.4297 0.5967 


b) 


Number 


of miRNAs 


1 3 5 8 10 15 30 


Average 


Number of FFLs 


1.51 4.56 7.65 12.32 15.45 23.21 46.48 


mi 


0.0186 0.0523 0.0920 0.1064 0.1356 0.1653 0.2949 


m 2 


0.0441 0.1138 0.1653 0.2052 0.2536 0.3111 0.4682 



Section [TTJ 

A. Derrida Values 

In this study, we introduce a new measure for the robustness of stochastic networks by 
quantifying the degree of determinism of network dynamics. Another measure that can be 
used in the Boolean context was suggested by Derrida in 1986 24 . Pairs of initial configura- 
tions of fixed Hamming distance are sampled from the entire state space, and their mean 
normalized Hamming distance, after being updated using update functions and propensity 
parameters, is denned as the Derrida value for a given initial Hamming distance. Lower 
Derrida values reflect more stable dynamics. To take time dependencies into account, we 
also considered the mean Hamming distance after two and three time steps as has been done 
earlier 25 . Table IVl displays the percent change in Derrida values starting with a basic 50-gene 
network and introducing 30 miRNAs. In all cases, the change is negative, i.e., the Derrida 
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TABLE IV. If each of a target gene's transcription factor regulates a miRNA that degrades the 
transcription factor mRNA, then only the fixed value a\ € {0, 1} can be taken on at a steady state 
because each transcription factor and its regulated miRNA have to take on the same value, or 1, 
at a steady state (gray rows). Here, the light gray rectangle contains the original update function, 
and ai, . . . , 04 are any Boolean values such that the update function depends on both inputs. 
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values decreased, indicating that the extended network exhibits more stable dynamics than 
the basic network, which we observed for different network sizes as well. Thus, another 
commonly used robustness measure also agrees with our hypothesis, which suggests that our 
findings are independent of the choice of robustness measure. 
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TABLE V. Derrida values for initial small disturbances of a Hamming distance up to 5 were 
simulated for a basic 50-gene network and an extended network with an additional 30 miRNAs 
of full strength. Multiple time steps were taken into account to consider time dependencies. The 
table shows that the percent change of the Derrida values from the basic to the extended network is 
always negative, indicating that our findings do not depend on the choice of the robustness measure. 



Hamming Distance 


1 


2 


3 


4 


5 


Time Steps 


1 


-2.8343 


-5.5589 


-7.7702 


-9.6424 


-11.1946 


2 


-3.6709 


-5.0620 


-6.3319 


-7.3764 


-8.3597 


3 


-6.4927 


-7.4783 


-8.3548 


-9.0794 


-9.8100 



V. Discussion 

We have examined the effect of feedforward loop motifs in stochastic Boolean network 
models of transcriptional networks, in analogy to the regulatory effects of miRNAs. Our goal 
was to test the hypothesis that these regulatory motifs have the effect of buffering the network 
against stochastic effects in the sense that they stabilize the basins of attraction. To this 
end, we conducted a computational experiment on a large number of randomly generated 
networks. The networks were modified by introducing additional nodes and feedforward 
motifs in a way that suggests regulation by miRNAs. To capture the effect on network 
stability we introduced a new measure of stochasticity of a network suitable for this purpose. 
Using this measure, as well as the classical measure of Derrida values, we showed that indeed 
the introduction of miRNAs has the hypothesized buffering effect. 

The number of miRNAs that are introduced strongly influences the magnitude of the 
stabilizing effect, so that one might wonder how many feedforward loops can be expected to 
be found in actual gene regulatory networks. In a data set from E. Coli, among 424 nodes 
with 519 edges, 40 FFLs have been found^S. In S. cerevisiae, among 685 genes with 1,052 
interactions, there are at least 70 FFLs 27 . However, restricting the data to subnetworks, 
we find other occurrence frequencies of FFLs. A subnetwork of E. Coli of 67 nodes with 
102 edges containing 42 FFLs was identified (some new FFLs had been found by then), 
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and in D. melanogaster, a subnetwork of 54 nodes and 167 edges contained as many as 157 
FFLs^. These numbers indicate that the question of how many FFLs are reasonable in 
a gene network of a certain size seems to depend strongly on the average in-degree of the 
nodes; whereas even large networks with average in-degree of less than 2 have few FFLs, this 
number can rise quickly when the network becomes more highly connected, as indicated by 
the considered network of D. melanogaster, with an average in-degree of approximately 3. 

Additionally, we looked at the correlations between the number of attractors in both 
networks and the number of common attractors, where we defined a configuration in both 
networks to equal if the states of all genes, i.e., the first n bits, coincide. Figure [2] shows 
the observed correlations, and we notice expected decreasing correlations between all three 
variables when more miRNAs are introduced. Surprisingly, the number of attractors of the 
extended network is much more strongly correlated with the number of common attractors 
than the respective number for the basic network, the cause of which remains to be explored. 

This study can be extended in several ways, which we are planning to pursue. To make 
the study design more realistic, it is useful to introduce additional nodes for proteins, in 
order to be able to implement more mechanistic details of miRNA regulation. Also, here we 
do not restrict the regulatory rules to those that correspond to activation and inhibition only, 
which does not allow the classification of feedforward loops into coherent and incoherent, 
an important distinction. Also, a more careful study remains to be done on the effect of 
miRNAs relative to their position in the network and the local network topology into which 
they are embedded. Finally, another limitation of this work is that only intrinsic noise is 
being considered as a perturbation. It is important, however, to also take extrinsic noise 
into account, which requires an extension of the SDDS framework. 

VI. Conclusions 

This study provides computational evidence that miRNA-mediated feedforward loops 
have the effect of buffering the network against phenotypic variation due to stochastic effects. 
Introducing a feedforward loop motif has a local effect on network dynamics that propagates 
to a generally much smaller global effect on attractor stability. Thus, as the number of 
feedforward loop motifs increases, the overall stabilizing effect increases as well. In our 
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study, the number of miRNAs introduced is of a relative order of magnitude that might be 
expected in an actual transcriptional network. Thus, our computational experimental setup 
can be used in conjunction with an appropriate experimental system to investigate the effects 
of individual miRNA actions. 
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FIG. 1. m\ and rri2 are plotted against the number of introduced miRNAs. Networks are not 
necessarily strongly connected, and the miRNA has full strength. The size of the considered 
networks varies from 5 (solid line) to 50 (dotted line). The impact of a single FFL on the dynamics 
is larger in smaller networks, which suggests that larger networks require more FFLs for the same 
amount of stabilization. 
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FIG. 2. The correlation between the number of attractors in the basic and the extended network 
and the number of common attractors is compared pairwise. These values have been obtained 
from not necessarily strongly connected networks of 50 genes and the introduced miRNAs had full 
strength. 
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